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Abstract. We show that for steady compressible potential flow in a class of 
straight divergent nozzles with arbitrary cross-section, if the flow is supersonic 
and spherically symmetric at the entry, and the given pressure (velocity) is ap- 
propriately large (small) and also spherically symmetric at the exit, then there 
exists uniquely one transonic shock in the nozzle. In addition, the shock-front 
and the supersonic flow ahead of it, as well as the subsonic flow behind of it, 
are all spherically symmetric. This is a global uniqueness result of free bound- 
ary problems of elliptic-hyperbolic mixed type equations. The proof depends 
on the maximum principles and judicious choices of comparison functions. 



1. Introduction and Main Results 

This paper is devoted to establishing uniqueness in the large for a class of tran- 
sonic potential flows with shocks in three-dimensional divergent nozzles. These 
transonic shocks are spherically symmetric and widely used in aerodynamics and 
computational fluid dynamics to simulate those transonic shocks appeared in the 
so called de Laval nozzles (i.e., convergent-divergent nozzles, cf. pi [51 \TU[ [TTl I15j ). 

To formulate the problem, it would be somewhat convenient to use (r, 9, (f>), the 
spherical coordinates of a point in R'^, with O<6'<7r,O<0< 2tt. Let S be a 
domain on the unit 2-sphere S^ C R^. For two fixed positive constants < r^, we 
call = {{r,d,(j)) I < r < r^,{9,(t)) e S} a straight divergent nozzle, and E its 
cross-section. The wall of the nozzle is the truncated cone T — {{r,d,(j))\r'^ < r < 
r^,{9,(l>) G as}. For i = 0, 1, = {(r, 6*, (/))|r = r\{9,4>) £ t} is called the entry 
and the exit of the nozzle respectively. Then obviously dft = F U S° U S^. Recall 
that the standard Euclidean metric of R'^ written in the local spherical coordinates 
is G = dr(E)dr + r'^{d9 (g) d9 + sin^ 9d(j) ® d(p). 

We consider steady potential polytropic gas flows in Q. The governing equations 
are the following conservation of mass and Bernoulli law (cf. [HIIH]): 

div(pgrad(y9) = 0, (1.1) 

1 p"^-^ - 1 

-(grad(p, grad.^) H — = bo- (1.2) 

z 7 " t 

Since we use spherical coordinates rather than Descartesian coordinates, here div, 
grad, (•,•) are the divergence operator, gradient operator and inner product with 
respect to the metric G respectively. The unknown ip is the velocity potential (that 
is, grad(/3 is the velocity of the flow), bo is the known Bernoulli constant determined 



Date: March 30, 2009. 

2000 Mathematics Subject Classification. 35J25,35B35,35B50,76N10,76H05. 
Key words and phrases, uniqueness, transonic shock, free boundary, Bernoulli condition, max- 
imum principle, potential flow, nozzle. 

1 



2 



LI LIU AND HAIRONG YUAN 



by the incoming flow and/or boundary conditions, p is the density, and 7 > 1 is 
the adiabatic exponent. (Similar to [6 , our results and proof also hold for the 
isothermal case 7=1.) The pressure of the flow p and the speed of sound c are 
determined by 



P{p) = -, c'{p) = p-^-\ (1.3) 
7 



Expressing p in terms of w = v^(grad(^, grad(^} ^ ^ {d^ipY + ^^j^ + j^r^^^- 

p^piv')^{l + h~l)ibo~lv^))^, 

equation (jl.ip becomes a second order equation for tp: 

div(p((gradv3, gradiys)) grad(p) — 0. (1.4) 

It is well known that this equation is of mixed elliptic-hyperbolic type in general 
[8j; it is elliptic if and only if the flow is subsonic, i.e., i; < c or equivalently, by the 

BernouUi law, u < c* := ^ (l + (7 — l)6o) . 

Remember that we call E*^, the entry and exit of the nozzle respectively. This 
means that, we always assume: 

(H) drf > on E% i = 0,1. 

That is, the gas flows in fl on and flows out of on E^. 

Tremendous experiments and numerical simulations have shown that, for a given 
supersonic flow near the entry of the nozzle, then by giving an appropriately large 
back pressure at the exit, a transonic shock must appear in the nozzle (see §147 in 
[5] or §4.3.4 in [H]): the flow is discontinuous across a shock-front; the flow ahead 
of the shock-front is supersonic, and behind of it is subsonic, with pressure increases 
(velocity decreases) across the shock- front. The position of the shock- front depends 
continuously upon the back pressure and vice versa. Therefore, to design nozzles 
work for special purposes in aerodynamics, one has to understand the existence, 
uniqueness, and stability of these transonic shocks via rigorous theoretical analysis. 

A basic strategy to attack this problem is as follows: (a) formulate a physically 
meaningful boundary value problem for the governing partial differential equations 
(PDE), and look for some special solutions involving transonic shocks; (b) study the 
stability of these special solutions under perturbations of the upcoming supersonic 
flow, the shape of the nozzle, and the back pressure etc.; (c) study uniqueness of 
the special solutions in a large class of functions. In step (a), instead of studying 
the full Euler system, one always constructs some simplified models (such as the 
potential flow equation, the popular quasi-one-dimensional model of nozzle flow, 
see [H [151 m] and references therein) , and then might obtain solutions with partic- 
ular symmetry by solving the reduced ordinary differential equations or algebraic 
equations. Step (b) essentially involves nonlinear small-perturbation problems: for 
example, one solves various linear PDE and then applying some elegant nonlinear 
iteration techniques to show existence of the nonlinear problem, see, for instance, 
[llElllHIIlIlOlllllIlllin]. In general, it is felt that step (c) is more harder, which 
involves nonlinear problems which are not of small-perturbation type. 

This paper devotes exactly to establishing the uniqueness in the large of a class 
of spherical symmetric transonic shocks in the divergent nozzle ft. Such spherical 
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transonic shocks were constructed in |8] and then widely used to explain the tran- 
sonic shock phenomena in de Laval nozzles (cf. [HITllIT^). Their stability, for the 
full steady Euler system, in two-dimensional case, was studied by '10] recently. See 
also [13, J4J . These results provide strong theoretical supports for the applications 
of these special transonic shocks in practices of aerodynamics, and in computational 
fluid dynamics to test various numerical schemes designed to capture shocks (cf. 

mum)- 

We remark that in [S] , Chen and Yuan have proved uniqueness of a class of flat 
transonic shocks in straight ducts, and as a byproduct, demonstrated that these 
flat shocks are instable, therefore not physical. This result is also consistent with 
previous instability results obtained in the papers [71 1121 116j , which were devoted 
to the stability issues by studying small-perturbation problems. 

Now let us formulate the boundary value problem of the potential flow equation 
in n and summarize some important properties of the special spherically symmetric 
transonic shocks. Then we will state our main result. The proof is left in the next 
section. 

We suppose that the flow is spherically symmetric and supersonic (i.e., v > c.^,) 
on E"; spherically symmetric and subsonic (i.e., v < c,) on S^. More specifically. 



e (c*, y 2(6o + 7^)) and a constant vi £ (0, c,), we consider the 
following problem: 

(fTID in n, (1.5) 

ip^u^r^, drLp = u° on S°, (1.6) 

(gradiy9, grad(/3) = on T} , (1.7) 

(grad^,7V} =0 on r, (1.8) 

where N is the outward unit normal of F. (If the outward unit normal of dY, is 
(ni, 712), then N = (0, ni/r^, n2/r^).) 

We remark that this problem is physically reasonable. In fact, since our purpose 
is to study transonic shock phenomena in 17, the flow should be supersonic at the 
entry, and the equation (|1.4p is hyperbolic there in the r-direction. So we need 
two initial- value conditions like (|1.6p . We choose (/? as in the first condition just to 
write the solution neatly. The flow is supposed to be subsonic at the exit, where 
(|1.4p is of elliptic type, hence one and only one boundary condition is necessary. 
Bernoulli type condition (|1.7p means that the speed of the flow is given at the exit, 
and by Bernoulli law, this is equivalent to impose a uniform back pressure at the 
exit, which is more physical than other conditions, such as itself (a Dirichlet 
condition) (cf. 8 ). In addition, (ll.8|l is the well known impenetrability or slip 
condition of inviscid flow along solid boundary. 

We will study uniqueness of solutions of (|1.5p - (|1.8p with the following structure 
(see [6]). 

Definition 1.1. For a function r = /(0, 0) defined on S, let 

n- = {{r,e,p)en\r<f{e,p)}, 

O+ = {(r,0,^)eO|r>/(0,^)}. 

Then (p e C°^^{^) n C^iyF) n C'^(yi+) is a transonic shock solution of (fTSll - lfrSl) 
if it is supersonic in fi^ and subsonic in il^, satisfies equation (jl.4p in il^ U 
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and the boundary conditions ()1.6|) - (ll.8|) point-wise, the Rankine-Hugoniot jump 
condition on S: 

p((grad(p+,grad(^+))(grad(^+,z/) = p((grad(^", grad<^"))(grad(p~, t/), (1.9) 

and the physical entropy condition on S: 

(grad(^+, grad(^+) < (gradiy9~, graders"), (1-10) 

where v is the normal vector of 5, and ip~^{ip~) is the limit value along S oi ip 
restricted in 51+ (il^). The surface S is also called the shock-front. 

For the existence of special transonic shock solutions and their important prop- 
erties, we have the following lemma. Without any ambiguity, we will also denote 
by the restriction of (p in 51+ (fi"). 

Lemma 1.1. Suppose the solution Lp depends only on r. Then for a given , we 
have the following results: 

(1) There is one solution [r) G C^(il) which is supersonic and solves problem 
(jl.Sp (|1.6p (jl.Sp in the whole nozzle f2. In addition, drip^ (r) is strictly monotoni- 
cally increasing on {r'^,r^). 

(2) There exists a connected open interval I C R+ such that for vi G /, there is 
uniquely one transonic shock solution 



/ \ I ^ V : ' ^ I' J ' s' 1 /I 1 1 ^ 

Mr) = < , ^ r 11 (1-11) 




to problem (fTSjl - fLSl) . with Sb := {(r, 6', </)) G 51 1 r = r^' e ir°,r^)} being the shock- 
front, (p^ (r) the supersonic flow obtained in (1), and f^{r) a subsonic flow. In 
addition, it has the following properties: 

(i) dr'p'^ir) is strictly monotonically decreasing on {rs',r^)- 

(ii) V3"(rs/) = (/3^(rs/) and (p~{r) > (p'^ (r) for r^/ < r < r^; 

(iii) drip^{r) > dr(p^{r) for rg' < r < r^; 

(iv) For any fixed R G (r*^, r^], drf^ (R) may be regarded as a continuous function 
ofrgi, and is strictly monotonically increasing for rs' G {r'^,R). 

Proof. Let v — dr^p and solve po from (-u°)^/2 -I- (pq~^ — l)/(7 ^ 1) = ^o- Consider 
the following algebraic equations of w(r) and p(r) for r G {r^,r^): 

r^pv = ao := {r^ p^u" , (1.12) 

'-v^ + P -i = 6o. (1.13) 

2 7 — 1 

The claims in the Lemma can then be shown by elementary calculus (see, for 
example, [3 [15]). In [15] the computation is carried out for the two-dimensional 
steady full Euler system, but the results there (i.e.. Propositions 1, 3 and Theorem 
6) still hold for three-dimensional potential flows (except the asymptotic expressions 
in Proposition 1 as r ^ cxd, which we do not need in this paper). Recall that v — dr(p 
here is the velocity of the flow in the r-direction. □ 

We also have the global uniqueness of supersonic flow Lp~ in the whole nozzle 51, 
which may be proved by standard energy estimates of nonlinear wave equations. 

Lemma 1.2. For a given , the solution •p>^(r) G C^(51) which solves problem 
p.5p p.6p p.8p is unique in the class of functions describing supersonic flows. 
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Now we state our main result. 

Theorem 1.1. Under the hypotheses (H), for a given u'^ and then any vi G /, 
there exists one and only one transonic shock solution to problem (|1.5p - (|1.8p in 
the sense of Definition [7771 which is exactly the (pbiT) constructed in Lemma \1.1\ 
corresponding to u'^ and vi . 

Remark 1.1. This uniqueness result also holds if we consider the whole spherical 
shell n' = {(r, P) |r° < r < r\P e S^} instead of the nozzle ft. The proof is 
similar and simpler. 

Since the supersonic flow ip~ is known, independent of downstream conditions 
and unique in the large according to Lemma 11.21 to prove Theorem 11.11 we just 
need to show that any possible shock-front must be Sb and the flow behind of it, 
(p^ , must be ip'^ . Therefore this is essentially a uniqueness result of a free boundary 
problem for a nonlinear second order PDE, with the transonic shock- front S being 
the free boundary. On the free boundary we have two boundary conditions, namely 
the Dirichlet condition 

ip+ = ip', onS (1.14) 
due to the continuity of tp across S (see Definition [TTT]), and the Neumann condition 

CH). 

In the rest of this paper. Section 2, we will present the proof of Theorem ll.il It 
might be a little unexpected to find that the proof is not very hard. It depends on 
maximum/comparison principles of elliptic equations and some judicious choices 
of the above constructed special solutions as comparison functions. The ideas are 
generalizations of those in ^ , and we need the deeper properties of the family of 
special solutions as in Lemma ll.ll to make them work. This fact, however, indicates 
that our methods might be applicable to a large class of free boundary problems 
which possess families of special solutions with fine structures. 

2. Proof of main result 

Let 

be a transonic shock solution to problem (|1.5p - (|1.8p . with S = {{r,9,(j)) G ri\r — 
f{9,4>) e C^(E), (61, (^) e S} being the shock-front. Then v, the normal of S, 
pointed from 17^ to Q'^ , in spherical coordinates, is 

(2.2) 

As in [6], by the Neumann condition (|1.8p and entropy condition (jl.lOp . we can 
easily show that 

(j^,iV) = ~-^(niae/ + n290/) =0, at 5 n F. (2.3) 

That is, the shock-front is always perpendicular to the wall of the nozzle. In addi- 
tion, direct computation yields 

p ■ (gradi^, v)^ p - (drp - ^depdef - ^ ^ d^pd^,/) . (2.4) 
V sm ' 
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We may write the equation ()1.4p in non-divergence form as 



2 



E 



^ A'^{DLp)d^ip + BiDip) = 0. (2.5) 



Here, for simpHcity, we set = r,y^ ~ 6, — <j>, and dij = dyiyj . We do not need 
to write out the specific expressions of (|2.5[) , but it is essential to note that A*-' and 
B depend only on the first order derivatives of ip. This can be seen from the well 
known form of (|2.5|) in Descartesian coordinates {x^,x'^,x^): 

3 

c^Aip - ^^^ip^^J^p^^^^Jip = 0. (2.6) 
«,i=i 

Indeed, we can also use this expression below to show the validity of maximum 
principles. 

Now suppose for the given vP and vi e /, the special transonic shock solution is 
tfib and its shock-front is {r = rs} n Cl. (See Lemma[TTTJ) The proof of Theorem ll.il 
is then divided into two cases. 

Case 1 . We first show that if there holds 

min/>r„ (2.7) 

s 

then f = Ts and w :— w{r, 9, (p) = ft^ir) — ^^{r, 6, (f) = 0. 

Let us consider the domain = {(r,0,(/)) e r2 1 r > fi6,(j))}. By (H^D, both 
and ip'^ are well defined in $1^. Therefore we may formulate a boundary value 
problem of w in as follows. 

First, w solves the following linear PDE: 

2 2 

i,j=0 i=0 
2 

2 

Y d.,p+{A^^{Dp+) - A^'{Dp+)) + B{D^+) - B{Dp+)) 
0. (2.8) 



-( 



Note that there is no zeroth-order term cw here. Since p^ is subsonic, this equation 
is uniformly elliptic, and by our requirements in Definition ll.il all the coefficients 
a'-', 6* are bounded. So by the strong maximum principle 9J, if w is not a constant, 

its minimum can only be achieved on the boundary d^l^ = U (fij n F) U iS. 
Second, by ()1.7p . we have the following boundary condition of w on E^: 

(grad(^^ -I- grad(^"'' , gradw) ~ v\ - vf = 0. (2.9) 

By assumption [H) and Lemma ll.li 

(gradv5+ -f grad(/7+, (1, 0, 0)} > 0. (2.10) 
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So (|2.9p is a linear oblique derivative condition to w. Similarly, on f2j n F, there is 
a Neumann condition 

(gradw, iV) = 0. (2.11) 

Therefore by Hopf boundary point lemma if w is not constant, the minimum of 
w also can not be achieved on U (fi^ n F). (For points on n {^~^ n F), as in 
[Sj, we may use a locally even reflection arguments to prove that a minimum can 
not be attained there, since the two surfaces meet there at a right angle.) 

Third, suppose the minimum is achieved at a point {R,P) € S (i.e., R = 
f{P),P e E). We note that w satisfies a Dirichlet condition here (cf. (|1.14p ): 

w = 5(0, 0) := ^+ifi9, 0)) - if {9, </.)) < 0. (2.12) 

The inequality holds due to property (ii) in Lemma ll.li with r^/ = r^. 

We first suppose that g admits a minimum at P, which is an interior point of 
E. Then dgg = d^g = at P and, since dr(p~{r) > dr'^'^{r) by property (iii) in 
Lemma ll.il (taking r^/ = r^), we also have 

defiP) = d^fiP) - 0, (2.13) 

which indicates that dgw{R, P) — d^w{R, P) — 0, and 

dgv+{R, P) = d^^+{R, P) = 0. (2.14) 

Indeed, by (|1.14|) we have (p'^ {f {9 , cj)) , 6 , (j)) — (p^ {f {9 , (j))) , hence, for example, we 
get {dr^+){def) + deip+ = drf^dgf, and therefore deip+{R, P) = by (jTT^ . 

Now consider the Neumann condition on 5. By (ITTg)) . ([^ . ((TTO)) and (PH)) . 
there should hold 

p(|9^(p+P)a^^+ = p{\drip~\^)drip~ (2.15) 

at (P, P). We may solve from this algebraic equation uniquely one dr^p^ < dr(p~ , 
which is not less than dr(p^ ■ (Note that R > Vg. Here we used property (iv) in 
Lemma ll.il with rg' — Ts and r^/ — R.) Hence we have drW — drf^ — drf^ < at 
(P, P) G S. However, we see that ly = (1, 0, 0) pointed into at (P, P) by ((TT3)) . 
this is a contradiction to the Hopf boundary point lemma, which asserts that there 
should hold drW > at (P, P), where w attains its minimum. 

For P e 9E, since 5* is perpendicular to F, then (|2.13p still holds, and the above 
analysis also works. 

Therefore w must be a constant in fl'^ . Now look at (|2.12p . we get dgf = d^f = 
for {9,(j)) £ E. So / is a constant. If / = r^, then w = and the uniqueness is 
proved. 

Now we show that / > Ts is impossible. Note that w is a constant and ip'^ 
depends only on r, so (p^ depends only on r for f < r < , and therefore 

f{r,0,(j3) ^\ +, . ... 1 (2-16) 
I tp^yr), J S f ^ T 

is a special solution to problem ()1. 5^ - 111. 8p . By Lemma TLII (iv). since / > Ts, we 
must have vi — drP^{r^) > drP^{r^) = fi, a contradiction as desired. 

Case 2. We now turn to the case that 

P, :=min/<r,. (2.17) 

E 
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We will prove by contradiction that this is impossible. 
Let 

^B{r) = < +f . p / / 1 (2.18) 

be a special transonic shock solution constructed in Lemma 11.11 for which {r — 
JJs} n is the shock- front. Set Vi = 9riy9^(r^). Then by property (iv) in Lemma 
11.11 (with i? = r^, Ts' = i?s and then rs' ~ rg), 

Vi < vi. (2.19) 

Consider the maximum of the function w :— (p'j^ — (p^ defined in fl'j . Now w 
satisfies a linear uniformly elliptic equation similar to (|2.8p , only with ip'^ replaced 
by (fi'j^. The oblique derivative condition on is now 

(grades + grad(^+, gradw) = -vj <0, (2.20) 

and the Neumann condition on n F is the same as ()2.1ip . So if w is not a 
constant, its maximum can only be attained on S, where the Dirichlet condition is 

w = g{9, 4>) ifi+ifiO, <P)) ~ V'ifi^, < 0. (2.21) 

By the definition of Rs, the maximum, 0, can be achieved at the point P G E where 
fiP) = Rs. Hence dgf{P) = d^f{P) = 0, as well as deip+{Rs,P) = d^^+{Rs,P) - 
and v{Rs,P) = (1,0,0), as shown in Case 1. 

By the Rankine-Hugoniot condition (jl.9p . drf'^ should satisfy 

p{\drLp+\^)drip+ = p{\drip^\^)drLp' (2.22) 

at {Rs, P). However, it follows from Lemma ITTT] that dr^^{Rs) is the only solution 
to this algebraic equation which satisfies the entropy condition. Therefore drW = 
at {Rs , P) , which is a contradiction to the Hopf boundary point lemma from which 
drW should be negative at {Rs, P) if w is not constant. 

Hence w is a constant and by p.2ip . it is zero. Now we obtain, from (|2.20p . a 
contradiction. Therefore the Case 2 is impossible. 

This finishes the proof of Theorem 11.11 

Remark 2.1. In the proof of Case 1, we can also obtain contradiction if w is not 
constant by analyzing the maximum as in Case 2. But in Case 2 we are restricted 
to considering only the maximum of w to obtain contradictions. 

Remark 2.2. The same proof also works for the uniqueness of cylindrical transonic 
shocks for potential flows in two-dimensional straight divergent nozzles. 

Remark 2.3. As noted in 6 , in the proof we just used the fact from Definition 1.1 
that the flow on the right hand side of the shock- front is subsonic (i.e., the entropy 
condition) . We do not need any assumption such as the flow should be subsonic in 
the whole domain S7^. 

Remark 2.4. The assumption {H) is only used to guarantee the condition (|2.10p 
at the exit to be oblique. Since for the family of special solutions ipi, constructed 
in Lemma ll. 11 dr(p'^{r^) is positive and bounded away from zero, say, dr<p^{r^) > 
Co > 0, therefore we can relax {H) a little by, for example, requiring only that 
dr(p > —Co/2 at the exit. 
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